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MATHEMATICAL  MODEL  FOR  DETERMINING  PATE  OF  KILL  OF 
MICROBIAL  AEROSOLS  BY  VAPOR-PHASE  DISINFECTANTS 


1.  INTRODUCTION 

The  purpose  of  this  research  was  to  develop  a  method  for  screening 
chemical  vapors  for  decontaminating  microbial  aerosols. 

The  first  task  was  to  develop  a  technique  for  determining  the  decay  rate 
of  microbial  aerosols  as  a  function  of  disinfectant  vapor  concentration  and  relative 
humidity  (RH)  •  To  accomplish  this,  a  simple  aerosol  chamber  system,  which  was  small 
enough  to  fit  in  an  exhaust-equipped  hood  and  which  could  be  conveniently  and 
reproduo Ibly  operated  by  one  person,  was  used. 

The  final  part  of  the  investigation  concerned  the  development  of  methods 
for  screening  disinfectant  vapors  by  applying  Raoult's  law.  Binary  solutions  of 
disinfectant  and  water  were  used  to  determine  the  maximum  decay  rate  of  a  microbial 
aerosol  in  a  disinfectant  vapor  at  a  given  RH.  Ternary  solutions  of  disinfectant, 
water,  and  a  nonvolatile  nonelectrolyte  were  used  to  obtain  any  selected 
concen t ra t i on  of  disinfectant  vapor  and  water  vapor.  The  relationships  between  the 
variables  were  reduced  to  equations  by  regression  analysis  and/or  described 
graphically . 

The  project  was  discontinued  before  its  coni  vision;  therefore,  most  of  the 
graphs  and  equations  are  based  on  minimal  data  and  some  relationships  were  not 
completely  defined. 

2.  MATER IALK  AND  METHODS 
2-1  Test  "  icti r La . 

Re  r  ra  t  fi  c.a  rvescens  (SM)  was  used  as  the  test  organism  because  it  is 
easily  cultured  and  identified  and  because  of  its  extensive  use  in  the  past  as  a 
simulant  for  vegetative  bacterial  agents.  Cultures  of  SM  were  transferred  weekly  by 
inoculating  50  ml  of  nutrient  broth  (Difco)  in  a  250-ml  Erlenmeyer  flask  with  a  few 
drops  from  a  previous  culture.  The  flask  was  incubated  at  27°  to  30°C  in  a  rotarv 
shaker-incubator  for  18  hours  and  then  stored  at  40. 

2 • 2  Small-scale  aerosol  chamber  technique . 

The  aerosol  technique  devised  for  this  study  consisted  of  aerosolizing  SM 
into  an  inverted  three-neck  600-ml  round  bottom  flask  (aerosol  flask)  and  then 
transferring  a  portion  of  this  aerosol  into  a  three-neck  5-6-liter  flask  (transfer 
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flask)  from  which  assays  of  the  viable  aerosol  population  were  made  periodically  to 
determine  the  exponential  decay  rate  according  to  the  formula: 

C  -  C0e~kt  (1) 

or 

InC  «  InC  -  kt 
o 

where 

C  Is  the  SM  aerosol  concentration  per  unit  volume  at  any  time  t 

t  Is  the  time  In  minutes 

Cq  Is  the  Intercept  or  C  at  zero  time 

k  is  the  proportionality  constant  that  will  be  referred  to  here 
as  the  decay  constant 

The  SM  aerosol  was  produced  In  the  aerosol  flask  by  12  squeezes  of  a 
rubber  bulb  connected  to  a  glass  Vaponefrln  nebulizer  containing  2  ml  of  SM  broth 
culture.  To  eliminate  any  concentration  effects  due  to  settling,  a  5-ml  syringe  was 
connected  to  the  nebulizer  reservoir  which  allowed  the  mixing  of  the  culture  In  the 
nebulizer  prior  to  aerosolization .  No  attempt  was  made  to  control  the  RH  in  the 
aerosol  flask  which  increased  rapidly  from  ambient  to  70%  to  90%  as  a  result  of 
aerosol  production  and  the  humidity  of  the  ambient  air.  The  aerosol  flask  was 
equipped  with  an  outlet  filter  to  allow  air  to  flow  in  or  out  to  atmospheric 
pressure  and  to  collect  the  residual  aerosol  after  completion  of  a  test. 

Two  minutes  after  aerosolization,  200  ml  of  aerosol  was  withdrawn  from  the 
aerosol  flask  with  a  syringe  and  transferred  immediately  into  the  transfer  flask. 
Measured  volumes  of  aerosol  (5  to  100  ml,  depending  on  the  expected  decay  rate)  were 
withdrawn  periodically  from  the  transfer  flask  with  a  syringe  and  immediately 
expelled  Into  the  intake  tube  of  a  340-hole  sieve  sampler1  drawing  room  air  at  its 
design  rate  of  28  liters  per  minute.  SM  particles  in  the  aerosol  samples  were 
Impinged  on  the  surface  of  tryptose  agar  (Difco)  plates  in  the  sampler.  The  plates 
were  incubated  at  30°C  for  at  least  24  hours  and  then  the  colonies  were  counted  and 
the  regression  equation  C  =  C^e-^  was  determined  by  the  method  of  least  squares- 
C  was  estimated  from  the  colony  counts  by  a  statistical  method  called  the  positive 
hole  method  by  Andersen.^  This  method  converts  the  colony  count  to  a  particle 
count;  e.g.,  339  colonies  =  339  positive  holes  «  2008  particles. 

Zero  time  was  taken  as  the  time  when  one-half  the  aerosol  volume  (50  ml) 
had  been  expelled  into  the  transfer  flask.  The  entire  100  ml  was  introduced  in 
0.1  minute.  In  syringe-sieve  sampling,  t  was  the  time  when  one-half  the  volume  in 
the  syringe  had  been  expelled  into  the  sieve  sampler.  The  time  required  to  expel 
the  entire  volume  into  the  sieve  sampler  ranged  from  0.02  minute  for  a  5-ml  syringe 
sample  to  0-1  minute  for  a  100-ml  syringe  sample.  The  sif ve  sampler  was  turned  on 
just  before  the  syringe  sample  was  introduced  and  turned  off  0.15  minute  after  the 
sample  was  expelled. 
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In  syringe-sieve  sampling,  the  minimum  value  for  t  was  0.30  minute  and  the 
minimum  interval  between  subsequent  samples  was  0.8  minute.  The  number  of  syringe- 
sieve  samples  per  test  ranged  from  eight  in  control  (no  disinfectant)  tests  to  one 
in  tests  with  decay  constants  k  greater  than  nine.  In  the  latter  tests,  the  extra 
data  point  needed  to  compute  k  was  the  intercept,  D  ,  from  control  tests,  nominally 
2  *  10^  SM  particles  per  milliliter.  The  maximum  volume  of  a  syringe-sieve  sample 
was  100  ml;  therefore,  the  limit  of  the  measurement  system  was  k  =  33  which 
represents  the  recovery  of  a  single  SM  particle  from  a  100-ml  sample  of  aerosol  at 
t  =  0.30  minute. 

No  dilution  accompanied  the  withdrawal  of  samples  from  the  transfer 
flask.  The  pressure  balancer  was  simply  a  rubber  glove  or  balloon  hanging  down 
inside  the  flask  which  expanded  into  the  flask  as  samples  were  withdrawn  from  the 
flask.  It  also  served  as  a  stirrer  although  experience  with  several  other  methods 
indicated  that  stirring  tends  to  increase  the  decay  rate  without  improving  the 
"goodness  of  fit"  of  the  regression  equation.  Given  the  size  and  shape  of  the  flask 
and  the  small  aerosol  particle  size,  adequate  mixing  was  probably  achieved  by  the 
slight  temperature  differences  on  the  sides  of  the  flask  and  the  fact  that  it  was 
moved  slightly  during  each  procedure  involving  it.  In  the  tests,  the  flask  was 
gently  swirled  momentarily  (0.05  minute)  before  each  sample  withdrawal. 

The  RH  of  the  air  in  the  transfer  flask  was  measured  with  narrow-range 
Aminco-Dunmore  lithium  chloride  electric  hygrometer  sensing  elements. 

2.2.1  Glass  syringes. 

Only  glass  syringes  greased  with  stopcock  grease  were  used  for  aerosol 
sampling.  Disposable  plastic  syringes  were  Cried  but  all  types  and  sizes  gave 
erratic  and  much  lower  recoveries  than  glass  syringes,  possibly  due  to  electrostatic 
charge  effects. 

2 . 3  Screening  methods  —  application  of  Raoul’-'s  law. 

2.3.1  Binary  solutions  of  disinfectant  and  water. 

The  saturation  concentration  of  a  disinfectant:  vapor  is  the  maximum  amount 
which  may  coexist  with  water  vapor  in  the  air.  Stated  another  way  —  the  saturation 
concentration  of  a  disinfectant  vapor  in  air  is  a  function  of  the  RH.  If  each 
component  obeyed  Raoult's  law,  the  disinfectant  vapor  saturation  concentration  would 
diminish  linearly  with  an  increase  in  the  RH. 


The  bactericidal  activity  of  a  disinfectant  vapor  at  saturation 
concentration  was  defined  by  determining  the  bactericidal  activity  of  the 
equilibrium  vapor  over  a  binary  solution  of  disinfectant  and  water.  Maximum 
bactericidal  activity  of  a  disinfectant  vapor  at  a  given  RH  would  be  expected  at  a 
disinfectant  vapor  concentration  at  or  near  saturation.  The  bacterial  aerosol  decay 
constant  k  at  disinfectant  vapor  saturation  concentration  will  be  referred  to  here 
as  =  maximum  value  of  k  at  a  given  RH. 


The  technique  consisted  of  sealing  50  ml  of  the  disinfectant  solution  in 
the  transfer  flask  containing  ambient  air  and  allowing  it  to  vaporize  for  at  least 
lb  hours  to  its  equilibrium  or  saturation  vapor  pressure  at  a  constant  temperature 
(nominally  25°C).  The  composition  of  the  equilibrium  vapor  is  a  function  of  the 
partial  pressures  of  the  components  in  the  vapor  and  the  total  vapor  pressure.  The 
partial  pressure  of  a  component  of  the  vapor  is  defined  as  its  mole  fraction  in  the 
vapor  times  the  total  vapor  pressure,  assuming  ideal  gas  behavior  (Dalton's  law). 

The  relationship  between  the  composition  of  the  vapor  and  that  of  the 
solution  is  described  by  Raoult's  law  which  states  that  for  an  ideal  solution  the 
partial  vapor  pressure  of  each  component  is  equal  to  its  mole  fraction  in  the 
solution  multiplied  by  the  vapor  pressure  of  the  pure  component;  i.e., 

p  =  xp  (2) 

whore 


p  =  partial  vapor  pressure 
p  =  vapor  pressure  of  the  pure  component 
x  =  mole  fraction  in  the  solution 

F.ven  for  a  nonideal  solution,  Raoult's  law  applies  to  a  component  as  its  mole 
fraction  approaches  unity  or  a  solution  may  show  ideal  behavior  over  a  limited 
concentration  range.  in  most  solutions,  there  is  some  type  of  influence  of  one 
component  on  the  other  so  that  equation  2  does  not  hold.  There  may  be  a  tendency 
for  chemical  action,  forming  a  compound  in  solution,  or  a  component  becomes  less 
volatile  in  the  presence  of  the  other  due  to  the  formation  of  nonvolatile  complexes 
or  there  may  be  dissociation  of  the  components  in  solution. 

For  a  nonideal  solution,  the  composition  of  the  vapor  can  be  estimated 
from  vapor-liquid  equilibrium  data  or  from  Raoult's  law  in  the  form: 

o 

p  =  Yxp 

where  y  Is  the  activity  coefficient  which  is  essentially  a  deviation  factor. 
Deviations  from  Raoult's  law  are  often  extreme  ns  the  concentration  of  a  component 
approaches  zero. 


Fxcept  for  lactic  acid,  the  disinfectant  vapor  saturation  concentrations 
were  not  assayed.  However,  estimates  of  the  disinfectant  vapor  saturation 
concent  rat  ions  would  be  useful  information. 

The  disinfectant  vapor  saturation  concentration  was  calculated  from  vapor 
pressure  data  given  in  the  literature  using  the  following  form  of  the  ideal  gas  law: 

X  =  M  xp°/RT  (3) 


10 


where 

X  =  disinfectant  vapor  saturation  concent ra t ion  in  grams  per  liter 
M  =  gram  molecular  weight 

R  *=  molar  gas  constant  for  ideal  gases  =  0.08205 

T  =  thermodynamic  temperature  in  degrees  Kelvin  (K)  =  273.15  +  °C 

If  the  activity  coefficient  y  is  known,  the  expression  becomes 

X  =  M  xp°Y/RT  (4) 

Vapor  pressure  p°  at  25°C  was  extrapolated  or  interpolated  from  a  compilation  by 
Daniel  R.  Stull'  or  from  the  International  Critical  Ashlcs4  using  the  approximate 
Clapeyron-Clausius  equation  which  states  that  a  pier  of  In  p°  against  T  1  gives  a 
straight  line.  For  example,  the  table  below,  from  the  International  Critical 
Tables,  shows  vapor  pressure  data  for  formic  acid,  one  of  the  chemicals  screened. 
Vapor  pressure,  p°.  is  given  in  nnllg . 


0 

p 

18.  9 

33.  I 

A  125.9  189.' 

2  7  C;  ,  e. 

198.1  552.1 

753.4 

760 

°c 

10 

20 

30  40 

8  0 

h0 

If 

w  0  90 

100 

100.75 

According  to 

the 

regression 

equat i on 

co  rip  u  ted 

f ron  t 

he  above  data, 

In  p° 

=  18.179 

4309.8  T  ‘i  the  vapor  pressure  ot  pure  anhydrous  ‘urmir  acid  at  25°C  (298. 15K)  is 
4  1  mmllg  =  0.054  atm. 

2.3.2  Ternary  solutions. 

Ternary  solutions  of  disinfectant,  waiter,  ’>-.d  a  nonvolatile  nonelectrolyte 
were  used  to  investigate  a  method  for  obtaining  any  selected  concentration  of 
disinfectant  vapor  and  water  vapor  by  mixing  the  appropriate  mole  fractions  in  the 
transfer  flask  and  allowing  evaporation  to  equilibrium  vapor  pressure  at  a  constant 
temperature  (nominally  25°C). 

3.  R KSUbTS  AND  DISCUSSION’ 

3.1  Sfl  aerosol  particle  size  distribution. 

Table  A— 1  shows  the  particle  size  distribution  of  a  control  (no 
disinfectant)  SM  aerosol  at  70%  RH  and  24°C  as  determined  with  a  six-stage  Andersen 
sampler.  Model  0101. ^  Contrary  to  the  recommended  procedure  previously  described, 
the  aerosol  in  this  particular  test  was  stirred  with  a  large  magnetic  stirring  bar 
rotating  at  approximately  100  rpm.  NMD  is  the  number  median  diameter  of  the  aerosol 
particle  distribution.  OSD  is  the  geometric  standard  deviation  of  the  distribution 
Indicating  the  degree  of  heterogeneity  of  the  aerosol.  A  perfectly  homogeneous 
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aerosol  has  a  CSD  of  1;  l.e.,  all  of  the  particles  are  the  same  size.  r2  is  the 
coefficient  of  determination  of  a  regression  line  indicating  the  "goodness  of  fit" 
of  the  data  points  about  the  regression  line.  A  perfect  fit  has  an  r2  of  1. 

SM  aerosol  particle  size  appeared  to  decrease  exponentially  with  tine. 
Regression  analysis  showed  that  the  predicted  NMD  at  zero  time  was  1.1  urn. 

Sieve  sampler  data  were  obtained  for  comparative  purposes.  Andersen 
sampler  recoveries  were  greater  than  those  with  the  sieve  sampler  because  the 
aerosol  RH  increases  rapidly  on  passage  through  the  Andersen  stages  favoring  viable 
collection.  With  the  Andersen  sampler,  the  SM  aerosol  decay  constant,  k,  was 
0.0556;  whereas,  that  with  the  sieve  sampler  was  0.0746. 

The  sieve  sampler  was  used  in  this  study  instead  of  the  Andersen  sampler 
because  it  requires  much  less  time  and  material  and  is  easier  to  use. 

3 •  2  Binary  solutions  of  disinfectant  and  water . 

The  results  with  several  binary  solutions  are  sunmmarized  in  table  A-2. 
The  data  are  limited  to  the  listed  RH's  for  each  rest  disinfectant,  thus  graphing  of 
kg  versus  RH  was  not  feasible. 

One  of  tile  formic  acid  solutions  tested  contained  90.6%  formic  acid  by 
weight.  Using  46.03  as  the  gram  molecular  weight  of  formic  acid,  the  composition  of 
tills  solution,  expressed  in  mole  fractions,  was  formic  acid  =  0.79  and  water  *= 
0.21.  Solving  equation  3,  which  assumes  ideal  behavior,  the  predicted  formic  acid 
vapor  saturation  concentration  would  be  0.080  gram  per  liter  at  25°c.  Similarly, 
the  predicted  RH  of  the  equilibrium  vapor  would  be  71"'.  since,  by  definition,  RH  = 
p/p°.  However,  formic  acid  and  water  are  known  to  form  an  azeotropic  mixture  with  a 
maximum  boiling  point.  Such  a  condition  is  caused  by  negative  deviations  from 
Paoult's  law;  l.e.,  both  components  have  partial  vapor  pressures  less  than  those 
predicted  by  Raonlt's  law.  Stated  another  way  —  the  activity  coefficients  y  are 
less  than  1.  Reduction  of  literature  datau,i  shows  that  the  equilibrium  partial 
vapor  pressures  of  formic  acid  and  water  over  an  aqueous  solution  containing  90.6% 
formic  acid  by  weight  are  0.039  and  0.0054  atm,  respectively,  at  25°C  which  converts 
to  an  activity  coefficient  of  0.91  for  formic  acid  and  0.82  for  water.  Therefore, 
the  predicted  formic  acid  vapor  saturation  concent ra t ion  according  to  equation  4 
would  be  0.073  gram  per  liter.  The  predicted  F.H  of  the  equilibrium  vapor  would  be 
0. 1 7  or  17%. 


Bactericidal  activity  of  formic,  acid  vapor  at  saturation  concentration  was 
much  greater  at  17%  RH  than  at  98%  RH.  However,  the  next  section  of  this  report, 
dealing  with  ternary  solutions,  shows  the  k  versus  RH  curve  for  formic  acid  vapor 
as  having  a  maximum  value  of  k^  at  about  75%  RH. 


Solutions  of  the  two  hydroxyca rboxyl Ic  acids  (lactic  acid  and  glycolic 
arid)  were  nonideal  in  that  they  contained  internolpcular  esters  in  equilibrium. 
The  lactic  acid  solution  equilibrated  was  an  aqueous  solution  87.9%  by  weight  total 


lactic  acid.  Lactic  acid  solutions  are  composed  of  lactic  acid,  intermolecular 
esters,  and  water  all  in  equilibrium.  The  esters  are  expressed  in  terms  of 
equivalent  quantities  of  lactic  acid.  Complete  esterification  would  correspond  to 
125%  lactic  acid  since  100  grams  of  a  polymer  (C^H^Oj)  hydrolysis  gives 

( 100)  (90)  /  (72)grams  =  125  grams  of  lactic  acid.  According  to  literature  data,*5  the 
approximate  composition  of  an  87.9%  by  weight  solution,  expressed  in  mole  fractions, 
is:  monomeric  lactic  acid  =  0.41,  intermolecular  esters  •=  0.08,  and  water  =  0.51. 

The  lactic  acid  vapor  saturation  concentration  was  measured  and  found  to 
be  3.7  x  10  5  grams  per  liter.  Solving  equation  3  and  using  90.08  as  the  molecular 
weight  of  lactic  acid,  the  vapor  pressure  p°  ot  pure  anhydrous  crystalline  lactic 
acid  would  be  2.5  *  10“^  atm  at  25°C.  Literature  data  on  the  vapor  pressure  of 
lactic  acid  are  sparse  and  in  poor  agreement  because  it  is  difficult  to  determine 
the  boiling  point  of  lactic  acid  at  reduced  pressure  and  impossible  at  atmospheric 
pressure  because  of  the  ease  with  which  it  undergoes  self-esterification  when 
heated.  Values  for  p°  at  25°C  extrapolated  from  the  literature  ranged  from  1  *  10 
to  5  x  10_&  atm. 

The  approximate  composition  of  the  glycolic  acid  solution  equilibrated, 
expressed  in  mole  fractions,  was:  monomeric  glycolic  acid  =  0.34,  intermolecular 
esters  =  0.04,  and  water  =  0.62. 

The  levulinic  acid  solution  was  actually  multicomponent  rather  than  binary 
In  that  the  solution  equilibrated  was  technical  levulinic  acid  containing 
lactones.  It  is  very  difficult  to  obtain  by  distillation  levulinic  acid  which  does 
not  contain  traces  of  lactones  and  water.  The  approximate  composition  of  the 
equilibrium  solution,  expressed  in  mole  fractions  was:  levulinic  acid  =  0.87, 

lactones  =  0.12,  and  water  =  0.01. 

Except  for  lactic  acid,  the  disinfectant  vapor  saturation  concentration  X 
was  calculated.  These  calculated  values  of  X,  particularly  those  where  the 
disinfectant  liquid  mole  fraction  x  does  not  approach  unity,  are  considered 
approximate.  According  to  Chemical  Engineers'  Handbook, 1  accurate  vapor-liquid 
equilibrium  data  are  difficult  tc  obtain  and  much  of  the  published  data  on  activity 
coefficients  are  of  doubtful  quality.  Chemical  Engineers'  Hardbook  also  describes 
several  graphical  and  mathematical  methods  based  on  thermodynamics  for  checking  the 
accuracy  of  activity  coefficients.  Published  data  on  vapor  pressures  would  be 
expected  to  be  more  reliable  than  those  for  activity  coefficients  but  should  not  be 
regarded  as  having  the  same  degrees  of  accuracy  as  that  for  the  most  widely 
researched  chemical,  water. 

3. 3  Ternary  solutions  of  formic  acid,  water,  and  glycerol. 

Ternary  solutions  were  tested  because  kg  was  too  large  at  most  PH's  to  be 
determined  with  binary  solutions  of  formic  acid  and  water;  i.e.,  dilution  with  a 
third  component  was  necessary  to  obtain  finite  values  of  k  at  most  RH's.  Also,  k 
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could  be  predicted  over  the  whole  vapor  composition  range  from  a  minimum  of 
experimental  data  If  a  functional  relationship  between  k  and  the  independent 
variable  (equilibrium  vapor)  exists. 

The  results  with  ternary  solutions,  assuming  Ideal  (Raoult's  law) 
solutions,  are  summarized  in  table  A-3.  Glycerol  was  used  as  the  third  component 
because  of  its  mutual  solubility,  lack  of  bactericidal  activity,  and  low  vapor 
pressure  (6. A  *  10  7  atm  at  25°C)  compared  to  formic  acid  and  water. 

The  quantities  listed  in  table  A-3  or  used  in  data  reduction  and 
regression  analyses  are: 

xj  “  formic  acid  solution  mole  fraction 
*»  water  solution  mole  fraction 

X3  =  glycerol  solution  mole  fraction  =  1  -  x,  ~  x0 

^2  =  formic  acid  vapor  density  (grams  per  liter!  assuming  an  ideal 

solution 


Xj/£X  =  formic  acid  vapor  density  fraction  =  Xj/Xj  +  X2  + 
y  =  vapor  mole  fraction  assuming  an  ideal  solution 

yj  =  formic  acid  vapor  mole  fraction 

0,0.  o  o 

-  x]P  j/xjp  j  +  x2p  2  +  x^p  3 

=  0.054x  i/0. 054x2  +  0.031x2  +  6.4  x  10“ ^3 

y®  =  formic  acid  vapor  saturation  mole  fraction 

=  formic  acid  vapor  mole  fraction  when  glycerol  solution  mole 
fraction  x3  =  0 

=  ( 1 -x  2  >  (0.054)/(1-x2)  (0.054)  +  0.031x„ 

1  ,  by, 

k  =  kGe  I 

kQ  =  k  for  control  (no  disinfectant)  SM  aerosols 

k  =  k  eb  y* 
s  o 

3.3.1  Curve  f ttt Ing. 

Linear  regression  analysis  (least  square  method)  of  the  data  at  43%  RH 
(table  A-3),  assuming  a  simple  exponential  relationship  between  the  variables  of  the 
form 


k  .  k0e»  n. 


showed  that 


k  at  43%  RH  =  0. 02 1 le 1 8 ’ 3y 1 


r  =  0.999 

ka  at  43%  RH  =  1.4  *  10J 


“T-* 


Figure  B-l  is  a  semilog  plot  of  k  at  43%  RH  versus  y^  showing  the  fitted  regression 
line  6.  It  is  not  possible  to  decide  whether  the  relationship  Is  statistical  or 
functional  from  looking  at  the  data  either  graphically  or  mathematically;  additional 
evidence  is  needed.  In  the  case  of  a  functional  relationship,  there  exists  an  exact 
mathematical  formula  relating  the  two  variables,  and  the  only  reason  that  the 
observations  do  not  fit  the  equation  exactly  is  because  of  errors  of  measurement  In 
the  observed  values  of  one  or  both  variables.  In  t.h  case  of  a  statistical 
relationship,  there  is  no  exact  mathematical  relationship  between  the  variables  - 
only  a  statistical  association.  The  distinction  is  important.  If  equation  3 
describes  a  functional  relationship,  then  k  should  be  determinable  over  the  whole 
vapor  composition  range  from  a  minimum  of  experimental  data.  Otherwise,  the  fitted 
equation  could  be  reliably  used  only  for  interpolation  and  to  summarize  the  data. 

Alternate  relations  for  decribing  the  aerosol  decay  rate  are: 

k  vs  x,  or  X,  (7) 

and 

k  vs  X,  /::  X  ) 

The  regression  equations  describing  relations  7  and  8,  assuming  a  simple  exponential 
relationship  between  the  variables,  are: 

k  at  43%  RH  =  0.C3 7 le 34 ' 8x i  (9) 

r2  =  0.984 

ks  at  43%  RH  =  0.0371e34’ P( ’ _X3) 

=  4.7  x  10s 
or 

ks  at  43%  RH  =  0.037  le348>'l  (10) 

since  the  density  of  formic  acid  vapor  at  25°C.  according  to  equation  3  is  0.10  grams 
per  liter  and 

1 2. ox  /:: x 

k  at  43%  RH  =  0.01 12e  1  (11) 

r “  =  0.984 
kg  at  43%  RH  =  161 

All  of  the  regression  lines  fit  the  data  welL.  However,  regression  line  6 
fits  the  data  points  better  than  either  regression  line  9  or  regression  line  11. 
The  value  of  kg  predicted  by  equation  9  appears  Inordinately  high.  The  value  of  k^ 
predicted  by  equation  11  is  less  than  any  observed  or  reported  value  for  control  SM 
aerosols . 

Figure  B-2  is  a  semilog  plot  of  k  at  43%  RH  versus  xp  showing  fitted 
regression  line  9.  Figure  B-3  is  a  semilog  plot  of  k  at  43%  RH  versus  Xj/r.X, 


showing  fitted  regression  line  11.  The  patterns  of  the  observed  values  indicate 
that  a  curved  line  with  a  decreasing  slope  would  fit  the  data  better  than  regression 
line  9;  whereas,  a  curved  line  with  a  positive  curvature  would  fit  the  data  better 
than  regression  line  11.  The  best  fitting  lines,  assuming  no  particular  functional 
or  exact  relationship  between  variables,  were  described  by  the  following  second 
degree  exponential  equations: 


k  at  43%  RH 

2 

r 

ks  at  43%  RH 


0. 02 37e52' 2x 1  ”  81,9xl 

0.998 

15 


(12) 


and 


k  at  43 Z  RH 

2 

r 


kg  at  43%  RH 


0.0237e 

0.998 

718 


5 . 5  3Xj /E  X  +  9.28(X  /EX) 


2 


(13) 


Figure  R-4  is  a  graph  showing  fitted  regression  line  13. 

Equation  12  must  be  an  approximation  useful  only  for  interpolation  because 
the  predicted  value  of  ks  is  15;  whereas,  the  observed  value  of  k  at  Xj  *  0-20  is 
31.  This  is  impossible  since  k  cannot  be  greater  at  x^  =  0.20  than  at  x^  =  0.47. 
Stated  another  way  —  differentiation  of  equation  12  shows  that  the  regression  line 
has  a  maximum  at  =  0.32;  whereas,  the  comparable  value  for  the  true  regression 
line  would  be  not  less  than  0.47.  Figure  B-5  is  a  graph  showing  fitted  regression 
line  12.  Figure  B-6  shows  additional  evidence  that  a  plot  of  In  k  at  a  constant  RH 
versus  X  has  a  negative  curvature  from  recent  work  vtth  a  single-phase  system 
(unpublished)  in  which  the  disinfectant  vapor  concentration  was  directly  measured  by 
experiment . 

3.3.2  Mathematical  model. 


It  was  concluded  that  the  mathematical  model  for  determining  the  rate  of 
kill  of  a  microbial  aerosol  by  a  disinfectant  vapor  at  a  given  RH  and  temperature  is 
equation  5;  namely, 


k  =  koeb  yi 


and  that  the  maximum  value  of  k  at  a  given  RH  and  temperature  is 


ks  =  vb  yi 


(14) 
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:>  contrast,  the  associations  between 


and 


k  and  X  j 
k  and  X^/2X 


are  appirenrly  statistical;  i.e.,  the  fitted  regression  lines  could  he  reliably  used 
■  o'.  f  >r  interpolat  ion.  Also,  it  is  unrealistic  to  conclude  that  a  relation  which 
describes  the  equilibrium  vapor  in  terms  of  weight  rather  than  number  of  molecules 
e on  !  d  he  f  line  t  i  ona  1  . 


Table  A-3  shows  that  kg  and  b  are  functions  ot  the  RH.  Figure  B-7  is  a 
semi  lop  graph  of  k  versus  RH  showing  a  maximum  at  about  /‘t.”  RH  at  25°C.  The  slope 
“efficient  b  of  regression  equation  5  increases  with  an  Increase  in  RH  as  shown  in 
figure  B-8.  Thus,  the  bactericidal  activity  of  lorn  1c  acid  vapor  at  25C'('  is 

■  '■'mp  1  e  t  e  ’  v  described  at  all  RH's  >43%. 


tiONCl.i'S  IONS 

The  small-scale  chamber  technique  would  he  useful  fur  screening  candidate 
disinfectant  vapors.  The  equipment  is  inexpensive  and  generally  available.  Tests 
'..dti  binary  and  ternary  solutions  are  easily  and  rapid  1  v  per  formed*  The  results 
would  be  expected  to  be  highly  reproducible  because  the  variance  associated  with  the 
independent  variable  (equilibrium  vapor)  would  be  inhorentlv  small.  It  is  possible 
to  predict  the  microbial  aerosol  decay  rate  over  the  whole  vapor  composition  range 
:r  or  limited  experimental  data  because  the  first-order  re  1  a t ionship  between  the 
rr’s.q  decay  constant  k  and  the  equilibrium  vapor  is  apparently  functional  when 
■  •  1 1  ■  •  r  comp  isit ion  is  expressed  thermodynamically;  i.e.,  In  mole  fractions.  Tn 
entrust,  the  relations  used  in  the  literature,  namely,  microbial  aerosol  decay  rate 
i  r  i  cc  nsr  ant  KH  .•■  rsrs  disinfectant  vapor  density  and  microbial  aerosol  decay  rate 
!'  I  constant  disinfectant  vapor  density  versus  RH  are  statistical  and,  therefore, 

.  etui  on  1  v  for  Int  orpelat  ion. 


R Kf OMNRNPA T I Oh’S 


The  in]  lowing  procedure  is  reeornnendec.  for  screening  candidate 

I  i  ■ :  n f  e o  t  ant  /a pors . 


.  !  A  candidate  vapor  would  first  be  tested  against  Ba c i 1 1 u  s  subtil  is  var 

nicer  spore  aerosols  of  known  particle  size  distribution.  Spore  aerosols  are  much 
n  re  resistant  to  decay  than  SM  aerosols.  A  vapor  effective  against  spore  aerosols 
would  he  expected  to  be  much  more  effective  against  other  microorganisms. 

..  The  initial  test  temperature  would  he  2S°r  or  other  constant  temperature 

. •  .  r  room  temperature. 


5.3  The  equilibrium  vapor  over  a  binary  solution  of  disinfectant  and  water 

would  be  used  to  determine  the  maximum  bactericidal  activity  of  a  disinfectant  vapor 

at  a  given  RH.  Tests  would  be  conducted  at  several  RH's;  e.g.,  0.02,  0.05,  0.10, 

0.30,  0.50,  0.70,  0.90,  0.95,  and  0.98.  These  data  would  be  reduced  to  a  kg  versus 

RH  curve.  Control  tests  would  also  be  conducted  over  a  wide  range  of  RH's.  These 

data  would  be  reduced  to  a  k  versus  RH  curve.  The  slope  coefficient  b  of  the  In  k 

o 

versus  y^  regression  equation  at  each  RH  would  be  determined  from  the  values  of  k^ 
and  kg.  These  data  would  be  reduced  to  a  b  versus  RH  curve.  These  curves 
completely  describe  the  bactericidal  activity  of  the  disinfectant  vapor  over  the 
whole  vapor  composition  range.  The  predicted  values  of  k  would  be  considered 
approximate  being  based  on  only  two  observations  at  a  given  RH;  i.e.,  kQ  and  kg . 
However,  the  information  would  be  adequate  for  screening  purposes  particularly  if 
the  curves  were  smoothed.  The  composition  of  an  equilibrium  vapor  should  be 
measured  directly  by  experiment.  It  is  possible,  by  the  application  of 
thermodynamics,  to  predict  vapor  composition  over  the  whole  range  of  concentration 
from  limited  vapor-liquid  equilibrium  data  and  also  to  check  the  accuracy  of  the 
predictions  by  applying  thermodynamic  checks  such  as  those  described  in  Chemical 
Engineers'  Handbook.  Also,  published  vapor-liquid  equilibrium  data  might  be  used  to 
predict  vapor  composition  if  the  thermodynamic  checks  indicate  sufficient  accuracy. 

5.4  Ternary  solutions  would  be  used  if  either  k  was  too  large  at  a  particular 
RH  to  be  determined  with  binary  solutions  or  more  observations  were  desired  for  the 
determination  of  the  slope  coefficient  b  of  the  In  k  versus  y^  regression  equation 
at  a  given  RH. 

5.5  Determination  of  bactericidal  activity  is  not  the  sole  criterion  of  a 
disinfectant  vapor.  The  questions  of  toxicity  to  humans,  carcinogenicity, 
corrosiveness,  flammability,  economy,  ease  of  use,  dissemination,  and  others  must 
also  be  answered. 
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*  NMD  =  number  median  diameter  in  urn;  CSD  =  geometric  standard  deviation 
**  r"-  =  coefficient  of  determination;  N  =  number  of  observations  used  to  compute 
regression  equation. 


Ambient  air  at  10%  to  20%  RH 


Geometric  mean. 

b  =  slope  coefficient  of  regression  equation 
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Figure  B-2-  SM  Aerosol  Decay  Constant  k  at  43%  PH  and  23  C  as  a  Function  of 
Fornlc  Acid  Solution  Mole  Fraction  x, 
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X.  IX 

Figure  B-4.  SM  Aerosol  Decay  Constant  k  at  43"  RH  and  JS  C  as  a  Function  of 
Formic  Acid  Vapor  Density  Fraction,  X,/;;x 
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Alexandria,  VA  22333 


Human  Engineering  Laboratory 
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Chief,  Bureau  of  Medicine  &  Surgery 
Department  of  the  Navy 
ATTN:  MED  3C33 
Washington,  DC  20372 

Commander 

Nava  I  Weapons  Center 
ATTN:  Technical  Library 

Code  343 

China  Lake,  pA  93555 

US  MARINE  CORPS 

Director,  Development  Center 
Marine  Corps  Development  and 
Education  Command 
ATTN:  Fire  Power  Division 

Quant  I co,  VA  22  134 

DEPARTMENT  OF  THE  AIR  FORCE 

HQ  Foreign  Technology  Division  (AFSC) 
ATTN:  TQTR 

Wr  I  ght-Patterson  AFB ,  OH  45433 
Comma  nder 

Aeronautical  Systems  Division 
ATTN:  ASD/AESD 

Wr  i  gh t-Pat t erson  AFB,  OH  45433 
HQ  AFLC/LOWMM 

Wr i gh t-Patterson  AFB,  OH  45433 

HQ,  A F SC/'SDNE 
Andrews  AFB,  MD  20334 

HQ  AMD/RC 

ATTN:  Chemical  Oefense  OPR 
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Air  Force  Aerospace  Medical  Research 
Laboratory 
ATTN:  AFAMRL/HE 
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Klrtland  AFB,  NM  87117 

1 

USAF  T  AWC/THLO 
Egl In  AFB,  FL  32542 

OUTSIDE  AGENCIES 

1  Battel le,  Columbus  Laboratories 
ATTN:  TACTEC 

505  K I ng  Avenue 
Columbus,  OH  43201 

Toxicology  Information  Center, 
WG  1008 

National  Research  Council 
1  2101  Constitution  Ave. ,  NW 

Washington,  DC  20418 

US  Public  Health  Service 
Center  for  Disease  Control 
ATTN:  Lewis  Webb,  Jr. 

1  Building  4,  Room  232 

Atlanta,  CA  30333 
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1  Stlmson  Library  (Documents) 

Academy  of  Health  Sciences 
Bldg.  2640 

1  Fort  Sam  Houston,  TX  78234 
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